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We consider a Bose gas with an attractive interaction in a symmetric double well potential. 
In the Hartree approximation, the ground state solution spontaneously breaks the symmetry of 
the trapping potential above certain value of the interaction strength. We demonstrate how the 
Landau- Ginzburg scheme of the second order phase transition emerges from the quantum model 
and show its link to the spontaneous symmetry breaking mentioned above. We identify the order 
parameter, the critical point and analyze quantum fluctuations around it. 



Ultra cold atomic gasses become ideal systems to in- 
vestigate many body physics both theoretically and ex- 
perimentally []| . Experimentalists can prepare quantum 
degenerate gasses in a wide range of trapping potentials 
and tune strength and even character of particle inter- 
actions [2|- Various phase transitions in ultra-cold Bose 
gases have been already considered 0, 0, Q • For example 
it is shown in Ref. [5| that a finite-size, translationally 
invariant, one-dimensional Bose gas with attractive in- 
teractions above critical value is fragile against a weak 
symmetry-breaking perturbation. The authors report a 
crossover from a translationally symmetric ground state 
to localized state. 

In this Letter we study a Bose gas with attractive in- 
teraction trapped in a symmetric double well potential. 
To the first approximation the ground state of such many 
body system is analyzed using Hartree approach, where 
it is assumed that all the atoms are in the same single- 
particle state (mode). This treatment leads to the non- 
linear Schrodinger equation, known as Gross-Pitaevskii 
Equation (GPE) @. The GPE is extremely successful in 
the analysis of systems of ultra-cold bosons and the sin- 
gle particle mode is called the "condensate mode" . When 
this equation is applied to the double well system, it turns 
out that the ground state solution of the GPE breaks the 
reflection symmetry of the double well problem, when the 
interaction strength exceeds some critical value. A better 
insight into the properties of the system is gained within 
the Bogoliubov method Q. In contrast to the GPE, it 
does not neglect the population of the non-condcnsatc 
modes, although their treatment is approximate. When 
the interaction strength in the double well system ap- 
proaches a critical value, the Bogoliubov theory reveals 
gapless spectrum, diverging number of particles in the 
non-condensate modes and simply breaks down. Hence 
in the vicinity of the critical value one has to analyze the 
properties of the ground state solution of the full quan- 
tum Hamiltonian. 

In present study we use tight binding a ppr oximat ion 
and apply the Bose-Hubbard model [^, M, 113, \lM- We 
show, using a simple approximation (called a continuum 
approximation [111 . Il2l |) that the Schrodinger equation 



obtained from our Bose-Hubbard Hamiltonian reduces 
to a one-dimensional Schrodingcr-likc equation of a ficti- 
tious particle in an effective potential. Evolution of the 
shape of this potential from a parabolic, through quartic, 
to double well as we change the interaction strength re- 
flects second order quantum phase transition. Finally we 
analyze quantum fluctuations in vicinity of the critical 
point. 

The Hamiltonian of atoms trapped in a double well 
potential reads 
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where oi (0,2) operator annihilates an atom in the left 
(right) well, Q stands for the tunneling rate between the 
wells and U is the on-site interaction strength. Here we 
consider attractive interactions hence U is negative. To- 
tal number of particles TV remains a conserved quantity. 
We can extract a constant part from the Hamiltonian ([T|) 
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which will be neglected in the further considerations. 

Let us review here some properties of the Hartree ap- 
proximation applied to the two-mode problem. It as- 
sumes all the atoms being in the same quantum state. 
The most general form of such a state is 
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Here, z is the relative population difference between the 
wells, (p is a relative phase and |vac) denotes the vacuum 
state. The Hartree method relies on calculating an ex- 
pectation value of the Hamiltonian ^ in the state ([3| , 
which gives 
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The dimcnsionlcss parameter 



7 = UN/n, 
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is a ratio of the on-site interaction per atom to the tmi- 
nehng rate. Notice that the minimum of the expectation 
value of the Hamiltonian ^ occurs for (^ = 0. In case of 
7 > — 1 the minimum appears for z = and for 7 < — 1 
it is shifted to the point izp, where zg = \/l — 7~^- The 
non-zero value of zq means, that solution of the Hartree 
approximation breaks the symmetry of the trapping po- 
tential. This suggests that a second order phase transi- 
tion takes place at 7 = — 1 in the corresponding quantum 
model. It is a main purpose of our investigation to con- 
firm this conjecture. 

As the total number of atoms is conserved, the wave 
function can be written in a Fock basis of states \N —n^ n) 
(here n denotes number of atoms in the right potential 
well), as 
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Then the stationary Schrodinger equation obtained from 
Hamiltonian ^ reads 

^ (- Un-HiA/(iV - n){n +1) + ^n-W{N - n + l)r 
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Let Zn denote the variable related to the relative popu- 
lation difference, 



_{N -n)-n _ 2n 

^" ~ N ^ " iV ' 

Using this variable, Eq.(l7]) transforms into 
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For a large number of particles the "■^" terms can be 
dropped [13| . Also in this limit the variable z„ may be 
treated as a continuous variable (i.e. Zn — > z) and we 
approximate the finite differences by a second order dif- 
ferential operation 
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Upon completing this procedure we reduce Eq. ^ to the 
following form 

^ ("- AyiT^^ + V{z)\ ^(z) = Ei;iz), (11) 



where we introduced the effective dimcnsionlcss potential 
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with 7 defined above. Notice that operator on the left 
hand side of Ea. (|lip is non-hermitian. However, our 
starting Eq.Q is a result of an action of a hermitian 
operator ^ on the wave-function ©. In principle, we 
could obtain the continuous limit in way conserving her- 
miticity, but the derivation would require much more 
sophisticated expansion and would loose its illustrative 
character. Instead we present an alternative derivation 
and we argue that the corrections required to restore her- 
miticity in Eg. pTj) are negligible in the large N limit. 

An alternative derivation employs the concept of rela- 
tive phase operator [ij]. This operator in general is ill- 
defined, but still is believed to give correct results in the 
large TV limit. We start from the classical Hamiltonian 
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where V{z) is defined by Ea. p2)) (see for example [15j). 
Next we perform canonical quantization by changing the 
conjugate variables z and ip into operators with the com- 
mutation relation [(p,z] = |^ [a, \l^. In the case of re- 
pulsive interactions, the phase fluctuations are dominant 
over fluctuations of population imbalance z. Thus in this 
case it is convenient to choose phase representation, i.e. 
replace z with ||--P- [ij|. On the other hand, in the case 

of attractive interactions the phase fluctuations are small 
cosi^ w 1 — ii^^ and it is more convenient to choose z 
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representation; (p ^ ^-r-. In this case we obtain 
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The Schrodinger equation takes the form 
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This is again a one dimensional Schrodinger-like equation 
for a particle in the potential U(z), where the coefficient 
2/A^ plays the role of h. One can compare Eq. (fTT|) and 
Eq. p^ and notice the the difference consists of two terms 
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To estimate the first term we notice that the characteris- 
tic width of the wavefunction ?/> scales like certain positive 



power of 1/7V (since 2/N plays the role of K). Hence each 
derivative over z generates a factor of the order of positive 
power of N , and so the first derivative is much smaller 
than the second one. The second term contributes to the 
potential and contains a factor of \/N'^. In conclusion 
both these terms can be neglected. 

Note that Eq. fTS)) turns out to generate accurate re- 
sults even for relatively low number of atoms. In order to 
evaluate the quality of the continuum approximation we 
compare the eigenstates (FigH]) and the energy level dia- 
grams (Figl2]) obtained by solving Eq. (|15[) and by numeri- 
cal diagonalization of the Hamiltonian ^ for 200 and 10^ 
particles. Already at the level of few hundred particles, 
the results obtained within the continuum approximation 
seems to be indiscernible from the exact ones. 

Another interesting issue is the relation between the 
expectation value of the full many body Hamiltonian in 
the Hartree approach and our continuum approximation. 
The potential V{z) in Ea. (jl2p is exactly the same as an 
expectation value of the Hamiltonian in Eq. ^ , provided 
we set (p = 0. The relation between them is established 
by the quantization procedure presented above. Recall 
that at the starting point of this procedure there was a 
classical Hamiltonian, which turns out to have exactly 
the same functional dependence as the expectation value 
of many body Hamiltonian in the Hartree approxima- 
tion (compare Eq.([3]) and Eg. p^ ). In the quantization 
procedure the terms with (1 — cos(y9) transform into the 
kinetic energy like operators, leaving V{z) as an effec- 
tive potential. Hence the same potential is responsible 
for spontaneous symmetry breaking in the Hartree ap- 
proach and for the quantum phase transition described 
by the many body quantum theory. 

In our opinion the real value of continuum approxima- 
tion is that one can predict the properties of the system 
by simply analyzing a form of the effective potential. Fig- 
ure [3] presents the shape of the potential for three crucial 
regions: before (7 > —1), at (7 = —1), and beyond 
(7 < —1) the critical point. In the first region the poten- 
tial has typically a quadratic form, while at the 7 = — 1 
it broadens substantially and becomes quartic. In the 
last region it has a form of double well with the minima 
at z = ztzg. This analysis paints a picture of the second 
order quantum phase transition, with the critical point 
at 7 = —1. 

The shape of the potential determines the shape of the 
ground state wave function, which is always symmetric. 
For 7 > —1 it is bell shaped and centered around z = 
with the width depending on 1/A^. For 7 < — 1 the wave 
function has a double hump structure, centered around 
±2o- This structure suggest the choice of an order param- 
eter. Our choice is the absolute value of the population 
difference between the wells: 
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FIG. 1: (color online). Comparison of the probability density 
[■(/'(zn)! obtained by diagonalization of the exact Hamiltonian 
([2]) (solid black lines) and by solution of Ea. (|15|l (dashed red 
lines) at the critical point, i.e. for 7 = —1. Panels (a) and 

(b) correspond to ground states while (c) and (d) to the ninth 
excited states. Results for A^ = 200 are presented in (a) and 

(c) and for iV = 10'' in (b) and (d). 




which in the continuum approximation is equal to A^|z|/2. 
In the classical theory the system would choose only 
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FIG. 2: (color online). Comparison of the energy levels (i.e. 
2(_E — Eo)/Q, where £"0 is the ground state energy) obtained 
by diagonalization of the exact Hamiltonian ([2} (solid black 
lines) and by solution of Ea. (|15p (dashed red lines) versus 7 
for iV = 200 (a) and iV = 10* (b). 



one of the wells. Hence it would be natural to define 
z as the order parameter. However in the quantum case 
the ground state of the system is the superposition of 
wavepackets localized in both wells. The modulus makes 
it insensitive to this superposition, it only examines the 
width of the wave-function. 

A clear sign of approaching the critical point is the in- 
crease of the fluctuations of the order parameter. The 
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FIG. 3: Effective potential (|12p for different values of 7 




variance of B in the ground state of the double well sys- 
tem versus 7 is shown in Fig. 3] for two different numbers 
of particles. Figure 2] indicates that the fluctuations are 
maximal around the critical point. With increasing value 
of N the width of variance decreases and its maximum 
tends to 7 = — 1. 



In conclusion, we analyzed the behavior of the ground 
state of the Bose gas with attractive interactions in a 
symmetric double well potential. Wc showed that this 
system experiences second order quantum phase transi- 
tion. This transition is revealed within the description 
using continuum approximation. This approximation al- 
lows for reduction of the many-body system in the two 
mode approximation (Bosc-Hubbard model) to a prob- 
lem of a fictitious quantum particle in an effective po- 
tential. Transformation of the effective potential from 
quadratic through quartic to double well with the in- 
crease of the interactions strength fits well to the Landau- 
Ginzburg scenario and allows to identify the order param- 
eter and the critical point. The presence of the critical 
point manifests itself by sudden increase of the fluctua- 
tions of the order parameter (so called critical fluctua- 
tions), which we confirm numerically. 



FIG. 4: (color online). Variance of the B operator, Ea. (|16|) . in 
the ground state of the double well problem obtained from ex- 
act diagonalization of the Hamiltonian ((2)1 (solid black lines) 
and by solving Ea. (|15p (dashed red lines) for A^ = 100 (a) and 
A^ = 1000 (b). 
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